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Abstract 



We develop a theory to calculate exciton binding energies of both two- and 
three-dimensional spin polarized exciton gases within a mean field approach. 
Our method allows the analysis of recent experiments showing the importance 
of the polarization and intensity of the excitation light on the exciton lumines- 
cence of GaAs quantum wells. We study the breaking of the spin degeneracy 
observed at high exciton density (5 lO^'^ cm^). Energy level splitting betwen 
spin -|-1 and spin -1 is shown to be due to many-body inter-excitonic exchange 
while the spin relaxation time is controlled by intra-exciton exchange. 
PACS number: 71.35 

I. INTRODUCTION 

The optical response of intrinsic semiconductors heterostructures has received consid- 
erable attention in recent years from both the theoretical and the experimental point of 
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view. The study of the luminiscence spectrum gives information about the lowest electronic 
excited state of a semiconductor, i.e. the exciton. Using polarized light, emited photons 
contain information on both exciton energies and their dependence on the exciton spin. |1| 
Time resolved photoluminescence experiments provide information on different exciton 
properties: exciton formation and decay processes, spin relaxation, binding energy evolu- 
tion, etc. In some of those experiments, performed in the picosecond range, |@-^ an energy 
splitting between spin +1 and spin -1 excitons has been reported. These studies show that 
the spin splitting increases with both the initial exciton density and the degree of initial 
polarization. 

So far, theoretical models have been proposed [|^,|[ to explain only the exciton spin 
relaxation without taking into account the exciton-exciton interaction. The models give 
alternative explanations to spin relaxation processes in terms of intraexcitonic exchange, 
PHT2| D'yakonov-Perel [jl3[ and Elliot- Yaffet [jl4[ mechanisms. However, these free-exciton 



models fail to describe the spin level splitting. 

As far as many-exciton effects are concerned, only the spinless high density exciton gas 
has been the subject of theoretical research in the last 30 years. |T5|-P^ Three different 
schemes have been developed. The first one, proposed by Keldysh and generalized by 
Comte and Nozieres, |jl6| is a BCS like approach. Another method, due to Hanamura and 



Haug, |jT5[ consists in writing the Hamiltonian in terms of exciton operators. When the center 
of mass momentum of each exciton is zero these two approaches are equivalent up to order 
(na'^)^ where n is the exciton densty, a is the exciton radius {fi)'^e/ fie"^, e being the electron 
charge, n the reduced exciton mass and e the dielectric constant, and d is the dimension of the 
space. The last availaible theoretical scheme consists in writing a Bethe-Salpeter equation 



and interpreting the homogeneous part as a multi-exciton Wannier equation. p!9| , |20[] The 



physics underlying the three approaches is always a mean field treatment of interaction 
between spinless excitons and so the equations obtained are analogous. The differences lie 
in the obtention of the equations and in the physical nature of the mathematical objetcs 
the theories are built up. In any case, spin splitting is beyond the scope of those spinless 



excitons theories. 

We present in this paper a theory of spin-dependent exciton-exciton interaction in two 
and three dimensions (2D and 3D). Such interaction produces in a gas with a difference in 
the spin populations, a level splitting. We concentrate on the spin and energy dynamics and 
their effect on the exciton recombination without paying attention to the exciton formation. 
Our theory for spin polarized systems is an extension of the exciton operator Hamiltonian 
approach (Hanamura-Haug approach [|l^). Our results for the 2D case are in agreement 
with experiments in quantum wells. Experimental work in the interacting regime with 
polarized light, remains to be done in bulk, to the best of our knowledge. 

This paper is organized as follows. In section II, we present the many-exciton Hamilto- 



nian 



\ T% from which we obtain many-exciton Wannier equations including spin. In section 
III we solve them using both perturbative and variational approaches. The main approxi- 
mation we make is to consider the exciton center of mass at rest. Our calculations, correct 
up to order (raa'^)^, give four exciton energies and wave functions as a function of density in 
2 and 3 dimensions. The maximum density experimentally reached is about 10^^ excitons 
cm~^. A reasonable estimate for the exciton radius is lOOA. When n = 10^^ cm~^ then 
na'^ is roughly 0.1. This means we neglect 10~^ compared with 10^^. Screening corrections 
to the energy within RPA are also taken into account altough they do not depend on the 
exciton spin. A rate equation is proposed and solved in section IV in order to obtain the 
time evolution of the different types of excitons. Once we have the energies as a function 
of densities and the densities as a function of time, we can get results directly comparable 
with experiments obtaining a qualitatively good agreement. Finally, section V is devoted to 
the discussion and summary of our main results. 

II. INTERACTING POLARIZED EXCITON GAS THEORY 

Dealing with exciton-exciton interaction is a complicated task. The exciton is a two-body 
excitation which is not, properly speaking, a boson. Even the non interacting exciton Green 
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function is not simple to calculate. pT|-p3| Consider the exciton creation operator: 



?/'| = dedh 0i(e, h)iplipl (1) 



where the index i is the set of exciton quantum numbers { K, J, M, z/}, K being the center of 
mass momentum, J the total angular momentum, M its third component and u labeling the 
internal state of the exciton. ipl creates a conduction band electron in e, where e = { r, Sg} , 
and analagously ip^ for a hole for which we neglect valence band mixing effects. On the 
other hand, 0j(e, h) is the exciton wave function. So, physically speaking, ipj is an operator 
that creates a electron-hole cloud following the excitonic probabilty amplitude. The exciton 
wave function can be written as a product of factors, namely, the center of mass term, the 
orbital part and the radial wave function. This factorization is possible if we assume that 
{ K, J, M, u} are good quantum numbers for the non-interacting exciton. 

Now it can be easily checked that exciton operators fulfill the following commutation 
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tpi^^pl = Si^i' - de dh (j)*{e,h){ dh'(f)i>{e, h')iljl^ljh + / de(f)i'{e' , h)iljlijje} (2) 



The second term of the left hand side of eq. ^ is an operator whose expectation value in 
a state with n excitons is of the order na'^. In the limit of vanishing density, the excitons 
are strictly bosons. Up to order {na'^Y it can be shown [T^ that the exciton operator 
Hamiltonian including the exciton exciton interaction can be written as: 



H = Hq + Hint (3) 



where 



Ho = Y.<t\T,+ n- V,hV > (4) 



and 



Hint = \ E 44i'i"i'^"'{<h^'\IdV\l"' > + <i,AhV\i"' >) (5) 



where 

< i|Te + Th- Veh\i' >= J dedh(j)*{e, h) {% + Th- Veh)(j)i'ie, h) 
< z, f >= J dedhde'dh'<f)*{e, h)<j)*,{e', h'){Vee' + - - ^he')0."(e, h)<f),,.{e', h') 

< I, z'mf, f >=-J dedhde'dh' 0*(e, h)<j)*,{e', h'){Vee' + - - ^he')0."(e', h)<f),,.{e, h') (6) 

It must be stressed that this Hamihonian is correct up to order {na'^y. This means that 
using this Hamiltonian we can (and we must) neglect the contributions to the energy of order 
{na'^y. The physical origin of Id is the direct unscreened Coulomb interaction between 
fermions belonging to different excitons while Ix is the inter-excitonic exchange, or the 
unscreened exchange interaction between fermions of the same type. It is important to 
distingish between inter-excitonic exchange and intra-excitonic exchange. The former is, 
as we stated before, a many exciton intraband (conduction-conduction or valence-valence) 
exchange and the later is a single exciton or interband (conduct ion- valence) effect. [p|-pn| The 
intra-excitonic exchange does not break the symmetry between spin +1 and spin -1 exciton 
and has a very weak influence in the exciton energy levels. |0 In this paper we neglect the 
intra-ecitonic exchange in the calculation of the exciton binding-energy. Nevertheless, the 
intra-excitonic exchange plays a very important role in the spin flip mechanism. 0,^ 

We use a mean field approximation. First, we calculate the expectation value of the 
Hamiltonian with a wavefunction equal to the product of the non-interacting exciton wave- 
functions: 

< H >=Y,^i + ^ Y,i<'i','i\^d + Ix\i',i> + <i',i\Id + Ix\i,i' >) ni> Hi (7) 

with 

00^ =<t\Te + n-Veh\t> (8) 

Now, if we make a functional derivation without any further assumption about the single- 
exciton wave function, the Euler-Lagrange equations we obtain are terribly complicated. 



15| , p!7|p^ One of the complications is non-locality, another is our lack of knowledge of rij. 



the quantum non equilibrium distributions. If we assume K = 0, say, all the excitons are 
at rest, nonlocality disappears and instead a set of continuous functions r;,(K, J, M, u) we 
have a set of discrete numbers n{J,M,u). This is the most drastic assumption we make. 
In the case of resonant excitation, the K = hypothesis is more realistic than in the 
non-resonant excitation case because, in the resonant regime, the system receives just the 
energy required to create the exciton without any kinetic energy. Experimental information 
is available both in the resonant and in the non resonant regime. We consider 

only the resonant case, i.e. K = 0. A less drastic and more usual assumption is that 
all the excitons are in their ground state u = 1. So, spin M is the only quantum number 
labelling our excitons. Since holes at the top of the valence band have an angular momentum 
3/2 and electrons a 1/2 one, excitons have angular momenta running from 2 to -2. Zero 
angular momentum does not play any role in the usual experimental configuration of light 
propagation along the growth axis in 2D systems. Hence, in order to describe the system 
we have only to know four numbers, i.e., the exciton populations which we denote by the 
array n = {n_|.2, n+i, n_2} = {n+2,n^,ni,n_2 }. Another approximation implicit 
to the derivation of the Euler-Lagrange equations is that for each n the energy takes its 
equilibrium value. This is equivalent to say that the changes in n are long compared with 
collision exciton times: the excitons interact between themselves many times before the 
populations change. This approximation is usually called adiabatic or quasi-equilibrium. 
Therefore, n can be taken as fixed in the theory to compute energy levels. 

Before obtaining the Euler-Lagrange equations, let us discuss the shape of the Hamilto- 
nian. Substituting z by M in eq. < H > is a sum of integrals, but all the terms < Id> are 
zero. This is a nice consequence of the approximation K = and the neutral charge of the 
exciton. We factorize the expectation values < > in spin part and spatial part. Let us 
first discuss the former and second the later. The spin part depends on the dimensionality 
of the space. We shall use the following notation: C,se,s,X^) probability amplitude 

of the electron having spin Sg and hole having spin Sh in an exciton with spin M. The spin 
part of < Ix > reads: 
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ImJ: = E ^ a,., (Ml) U,sAMs) Us'SM.) (9) 

Obviously the spin wave functions depend on the confinement of the excitons. We shall 
make the following approximation for the 2D case 

For the bulk case ^se,Sh(M) is set equal to the Clebsch-Gordan coefficient with J = 2 in 
analogy with the 2D case in which third components of the angular momentum equal to 
±2 appear. If the actual 3D exciton wavefunction should be a combination of J = 2 and 
J = 1, some quantitative changes would appear although the final results would remain 



qualitatively unaltered. Using the eq. |T0|, the only non-zero spin terms for the 2D case are: 



^2D = II'} 

Hi = it-i = = m = 1 (11) 



as well as the ones generated by the following symmetry properties: 



(12) 



In the 3D case the only non zero spin terms are: 

7-2,2 _ 1 4 _ 7-1,1 _ 10 7-1,-1 _ 6 
-'2,2 — -L ^3D = ^1,1 — Ye -'1,-1 ~ 16 

7-2,1 _ 3 7-1,2 _ 1 r2,-l _ 1 

-'2,1 ~ 4 -'2,1 — 4 -'2,-1 ~ 4 

and the ones generated by eq. [T^. The main difference between equations |TT] and |T3] is that 



(13) 



fd = IZl'J is not zero in the 3D case. These spin term are proportional to the interaction 
term ^m^-a/ in eq. |^. The physical consequence of this fact is that, in 2D, M-excitons do 
not interact with (— M)-excitons and this enhances the splitting when n{M) ^ n{—M). As 
it will become clear in section III, if we had IZl'l = ll'l then the splitting would be zero. 
In order to compare with the experimental results performed in quantum wells, we shall 



use our 2D results including a renormalized "/d term" (/Zi'i) as well as a non-interacting 
exciton binding energy E^^ = = 2Eq^. [^,0 We take fqw — 6/160, i. e. one 

tenth the 3D value. 

In order to work out the spatial part, we proceed to make the functional derivation: 
{6 < H > /6(f)M{r))/n{M) = to get four Euler-Lagrange equations, one for each M. 
In these four equations there are two well differentiated parts, one corresponding to 6 < 
Hq > /S(t>M{i") and the other corresponding to 6 < Hint > /S(t>M{i"). The former generates 
the usual Wannier equation [|l7l and the later generates the interaction terms. In order 
to reduce the four interaction terms V{e, e'), V{h', h), V{e, h'), V{e', h) to one, i.e. V(q), it 
is convenient to work in the momentum representation. We are going first to derive the 
equations in the case n(2)=n(-2)=0. This will simplify considerably the equations and will 
shed some light over the underlying physics. Besides, in the experiments the ±2 optically 
inactive |^ excitons are less populated than the optically active ±1. Hence, we have only 
two kind of excitons, namely, up and down, ] and [. We omit the algebra because it does not 
give any physical information. The two multiexciton polarized Wannier equations are: (we 
set ^ = 1) 

,2 



<r_ 

2/i 



Ef 



0t(9)-7t(9)+ 2Adn| |0t(?)I 7t(9) - Stt(9)0t(9) 



+ 







(14) 



2/i 



^1 



- ii{q) + |0i(g)r7i(g) - ^^^^(g) 



+ 



2/d n| 7T(9)0t(9)0i(g) + \h{<i)?ll -^]A<l)(t)l{<l) -^M{q)h{<i) 



(15) 



with 



f dt 

ls{q) = j 7^0s(r+gmi) 



(2vr) 

<P:{t + q)Mt + q) V{t) 



(16) 



where V{s) is the Fourier-transform of the bare Coulomb potential and E^ and Ey are the 
Lagrange multipliers associated to the normalization constraint which are interpreted as the 
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exciton energies. Note that equations 0, |T5| and depend on the dimension through 0, 
and fd- The first two terms in equations |1^ and |15| are the usual Wannier interactionless, 
two band, exciton equation. The other terms proportional to and n-^ represent the mean 
field exciton-exciton interaction. Observe that we can obtain each equation by reversing all 
the spins of the other because there is no magnetic field. The spin symmetry breaking comes 
from the fact that um 7^ n-u- 



III. CALCULATION OF THE INTERACTING EXCITON ENERGY LEVELS 



A. Perturbation theory 



Solving the multi-exciton Wannier equations 14 and exactly is not possible and there- 



fore we use some approximations. First we calculate E perturbatively, using the the two first 



terms in equations |TJ and |T5| as an unperturbed Hamiltonian and the interacting terms as 
the perturbation. The perturbation parameter is obviously na'^ and, as we stressed before, 
we must drop all the contributions to the enegy with order higher than na'^. Consequently, 
we do not go further than first order perturbation theory which means that the interacting 
exciton wave function is the same that the non-interacting one. Hence, to first order in 
perturbation theory, (j)i = (j)^ = where: 



(27r)5a 



[1 + (ag/2) 



21 1 



3^ . ofyra^^a 



are the exciton wave functions in the isotropic, parabolic two band model, the momen 



tum representation the Schrodinger equation is an integral equation. [28] The perturbation 
gives a first order correction 



Using eq. |lj in eq. |Tg we obtain the following expression for the 2D case: 
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where 



A^t = 2n| (Ji -12) + 2 fQw niih - h) 
AE^ = 2n^ ih -h) + 2 fgw n^ih - h) 



and 



Since we are in the lowest order of perturbation theory Ji and I2 do not depend on 
the perturbed wavefunctions. In the appendix we show that in 2D Ji = 7ra^|£'Q^| and 
J2 = vr|£;2i?|a2 3i5ji^ \El^\ is the 2D Rydberg El^ = -2eVa = -2h^ / ^a'^ = AE^^. 

We obtain 

AEt = /c|E2^| a' (n^ + fQw ni) 
AE^ = k\E^,^\ (n^ + fQw n^) 
A^^ ^ AE^ - AE^ = k\E',^\{l - fQw){n^ - n^)a' (20) 

where k=3.03. We observe that the inter-excitonic exchange interaction produces a blue- 
shift of the levels. Although this calculation does not include screening effects, A^^ gives 
properly the spin splitting because screening is spin independent as discussed below. We 
observe that A"^^ is proportional to the polarization (P = (n^ — ni)/{n-i + '^i)) and takes 
his maximum value when fqw = 0, i.e. in the strictly 2D case. Anyway, fgw is almost 
negligible compared to 1 and we shall drop it in the following calculations. In the non 



polarized case we retrieve the result obtained by Schmitt-Rink et al. for energy shifts. p9 
In 3D we obtain analogously Ii — I2 = 137r/3 which brings to 

AE| = 137r/3 \E^^\ + 



,16 ' 16 

AE^ = 137r/3 l^o'^l a' + ^ 

^ AE^ - AE^ = 13n/3\E',^\ Q - ^) (n^ - n{)a' = 3.4 \E',^\in^ - n^)a' (21) 
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The energy splitting between spin —1 and spin +1 excitons that occurs in the 3D case is, 
at equal densities na'^, very close to the 2D one both of them scaled in their corresponding 
Rydbergs E^. Hence, we predict that the energy splitting, measured in units of , has 
a very weak dependence on the quantum well width. On the other hand, in the 3D non 
polarized case we do not recover exactly the result obtained by Haug and Schmitt-Rink W7\ 
for energy shifts because we do not use spinless wave functions as they do. The effect of the 
spin part of the wave function is to set A^d = 10/16 instead of 1 as in their calculations. 

We have been able to calculate the many-body corrections to the exciton binding energy 
in the case n+2=n-2 = using perturbation theory. The many body corrections in E^ and 
E^ do not depend on the energy levels of ±2 excitons. Therefore, we do not need to write 
the multi-exciton equations for the n+2 7^ ^-2 7^ case. In order to avoid tedious algebra 
we can evaluate the energy corrections by counting how many integrals are non zero in the 
sums of eq. |^ noticing that each bracket makes a contribution equal to Ji — I2 times the 
spin factor. In this way we obtain for the 2D case the following interexcitonic exchange 
corrections 



-C/+2 



5 



E-i 

E-2 



I 



^0 « 



^+2 



(22) 



k k k ^ 
k k k 
k k k 
k k k 

A prediction of this theory is that, neglecting /Qvy-like factors, the M-excitons interact with 
equal strenght with all the others except with the (-M)-excitons for which no interaction 
exists. The energy splitting of ±1 excitons is given by 



V 



n+1 



A 



2D 



k\EQ |a (n+i - + (n+2 - ^^-2)- 



(23) 



The ±2 exciton density is negligible compared with that of the ±1 excitons because the ±2 
excitons cannot be optically generated in one photon processes. Therefore, their influence 
on A^^ is very small. 

In the 3D case we have 
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^ F ^ 



137r 



1 10 6 1 



-c/n a 



16 16 4 
4 16 16 

i 1 1 



1 



(24) 



This equation does not carry new physics compared to eq. ^ and, as in the 2D case, the 
influence of the ±2 excitons is hmited because of their very low occupation. In a situation 
with ri-i-i > 10^°cm~^ and ±2 exciton neghgibly populated, the later energy levels are higher 
(less bound) than the ±1 excitons both in 2D and 3D because the ±2 excitons interact with 
two higly populated excitons, say, +1 and -1. 



B. Screenig corrections 

Hitherto we have used the unscreened Coulomb potential reduced by the dielectric con- 
stant e of the material in its ground state. The presence of a considerable amount of excited 
mobile carriers (electrons and holes) screens the interacion between these carriers and the 
rest of the lattice. The renormalization of the exciton binding energy caused by screenig has 



received considerable attention. |[16|j20|j30|| In order to simplify the theory we have calculated 
the screening corrections in the Random Phase Approximation (RPA). In this approxi- 
mation the screening correction to the binding energy does not depend on the exciton spins 
T3,|D| being 



i\ <r ^ e|p«'^9 r _ „-i(3q r\ II \\2 1^1 olp-jag r _ i/3q r\ I ^ |\2 

AE- ^ -nj: V\q) ^-^^^ ''^^'^ (25) 

where M^x = + f^^h is the total exciton mass and | v > are the exciton internal states, 
a = ruh/MEx and (3 = rrie/MEx- Using the completeness relation we get: 

(< ls|l — cos{q r)\ls >)^ 



{Eo + 



2Mex' 



-2n Y.V\q) 



1- ' 



1 + a? q'/Af 



2 



^ (26) 



^0+ 2M^ 
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Next, we transform the summation in an integral following the usual prescription. The 
screening correction we obtain in the 2D case is: 



TT na'\E^''\w'F{w) 



(27) 



where w = (M/4/i)2 and 

F{w) 



°° dx 

X 



1 - 



(1 + a;2)3/2 



x^ + 



In a GaAs quantum well z = — ^ 2.5 which leads to w"^ 



42 



(28) 



1.2, w^F{w) ~ 0.41 



and AE^^ 



0.4l7r na'^\EQ^\. The screening contribution reduces the effect of the bare 



Coulomb interaction between the carriers producing a relative red-shift of all the exciton 
levels and therefore does not contribute to the spin level splitting. It must be stressed that 
w'^F{w) is a very smoth funcion of mh/me and the screening correction is rather insensitive 
to variations of this mass ratio. 

A physical interpretation of the way the screening modifies the exciton binding energy 
is the following. In the previous section we calculated the exciton binding energies taking 
into account the inter-excitonic interaction with the bare Coulomb potential Vq. Now we 
calculate the dressed Coulomb potential Vs in the RPA and we treat the difference Vg — Vq 
as a perturbation. It happens that, to the lowest order, Vs — Vq is proportional [|l^ to n 
and we must use the non-interacting wave functions to evaluate the screening correction. 
Furthermore, we cannot calculate the screening corrections to the blue-shifts caused by the 
inter-excitonic exchange interaction because they would be, at least, of order (raa^)^. Hence 
the 2D exciton levels including inter-excitonic exchange and screening are obtained from 

/ 



\E, 



2D I 



-1 
-1 



k — q k — q k — q —q 
k — q k — q —q k — q 
k — q —q k — q k — q 
—q k — q k — q k — q 



\ ( \ 



(29) 



where k = 3.03 is the bare interexcitonic blue-shift constant and q = 1.28 is the screening 
red-shift constant, both dimensionless. 
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The screening correction energy in bulk was calculated in RPA by Zimmermann ||30 
obtaining: 

. ^.r irnSD, f w(32 + Q3w + AAw^ + llw^ 8wU + 3w)\ o , , , 

AE- ^ ME'oV [-^ (iTW^) " /(^) (30) 

As in the 2D case, the screening correction in the RPA is quite insensitive to variations 
of z. Hence, we also take z = mh/nie — 2.5 for 3D GaAs wich leads f[w] = —5.0 and 



AEq'^ ~ — 57r|i?Q^|na'^. Hence, the bulk interacting exciton levels are given by eq. ^ minus 



the screening correction AE^ 



sc 




C. Variational Approach 



We have a set of complicated equations and |T5D for which we have applied first order 



perturbation theory. We cannot go beyond first order due to our previous hypothesis but we 
would like to extract more information from those equations. In order to do that we have 
tried a simple variational approach in the 2D case. As we did in the perturbation approach 
we treat eqs. |1^ and |T5| as a Schrodinger equation. We identify a Hamiltonian and minimize 

< (f){q,a)\H\(f){q,a) >, with 

0r(g) = 3 , (31) 

[l + (aag/2)2]f ' ^ ^ 

i. e. we use the exciton radius as a variational parameter. This Ansatz may be improved if 
we make a dependent on the spin. However, the calculations are much simpler with Ansatz 



31| . The non trivial integrals we have to perform are precisely Ji and I2 scaling a with the 
variational parameter a. After some algebra we arrive at: 

1 2 



EMia,n,nM) = \Eq^\ 



— h 3.03 UMO^a — l.28na'^(aY 



(32) 



where n is the total density and um is the M-exciton density. In this expression we have 
set /=0. Now we have to look for the a that minimizes EM{c(,n,nM) for each n. We have 
done that numerically obtaining that, up to na^ ~ 0.2, the variational technique and the 
perturbation theory predict the same energy AE(^M, n) with an error less than 1% and a 
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does not differ from 1 (the perturbation theory value) more than a few percent. We can 
conclude that, in the small density limit [na"^ < 0.2), the energy is properly given by first 
order perturbation theory and the wave function is the independent exciton one. 



IV. EXCITON SPIN DYNAMICS 



One could expect our theory to predict new spin flip channels originated by the in- 
teraction. However, this is not the case because the interaction terms in eqs. |l^ and |15 
are proportional to n and the transitions rates are proportional to the squared interactions 
terms, i.e. n^. Hence, following the considerations made in section II, we mfisi neglect these 
'interacting' transition rates: our theory predicts no significant variations of transitions rates 
with respect to those of the non-interacting theory while the energy levels are the corre- 
sponding to equation B^. Therefore, we borrow the population evolution equation from ref. 



— = W n 

dt 



where 



w 



-1,1 



V 







ex 







2,1 



wt 



2,-1 



(33) 



. (34) 



This equation is obtained from the master equation by making the approximation of quasi- 
equilibrium and taking into account the detailed balance principle. There are three 
kinds of transition rates: the easiest to understand is the radiative recombination rates Wr 
which affect only to the optically active excitons ±1. Also, we have the Wlj%t rates, say, 
the transition (M)-exciton to (M')-exciton caused by the e(h) spin flip. The last type of 



15 



transition rate, W^'^' is associated to the intra-excitonic exchange mechanism. Following 
ref. [0 we set 



W 



e{h) 
M,M' 



Wh = 1/AOOps 
1 



r^(^) 1 + exp[{EM' - E 



' ' ex 



1 



(35) 



l + exp[{EM' - Em)(3] 

where t'^''^ are the single particle spin flip times, [|7|,|8|,pr3| r'^^ is the exchange spin flip time 



calculated by Maialle et al., /3 = 1/KbT and Em are those of eq. The numerical 
values of r^^'^-' and r^^ are taken from the case II of the reference. |^ 



V. RESULTS AND CONCLUSIONS 

The solution of the non-linear equations ^ is obtained numerically by a Runge-Kutta 
method. Inserting the time dependence of the densities in eq. ^ we obtain the time 
evolution of the exciton levels. In our theoretical calculations the natural energy scale is the 
Rydberg (the 2D or 3D ). In order to compare our results with experiments ||2|-§ we shall 
plot the energies in units of Eq^. On the other hand, we have checked that the results are 
insensitive to the variation of the initial populations of ±2-excitons provided they are less 
than 10% of the total initial density , no and, consequently, we present flgures obtained from 
an initial density of ±2-excitons equal to zero. We adopt the following conventions: i) the 
(+l)-exciton is the more populated state at t = 0; ii) the populations are measured in uq 
units; iii) the origin of energies is taken at the bottom of the conduction band. 

In flgure |I] we plot our theoretical predictions of populations (a) and energy levels (b) 
with an initial density of na"^ = 0.1 (about 10^^ excitons per cm^) and a initial polarization 
P = (^1 — ni)/{n-^ + rii) = 80%. In flgure |l|(a) we observe that the polarization disappears 
in roughly 50ps. As we see in the inset the ±2 populations are, at least, one order of 
magnitude smaller and they decay much slower than the optically active ones. In flgure 
|l](b) the more important features are the spliting between spin +1 and spin —1 excitons and 
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the fact that the ±2-excitons levels are closer to the conduction band than the ±1 ones. 
In t=0 the splitting takes its maximum value, 0A5E^^ , about 3 meV, and then decreases 
becoming zero at t=50 ps. It must be stressed that the copolarized (more populated) exciton 
level shows an energy blue-shift and the contrapolarized (less populated), when is very low 
populated, starts having a red-shift and, when its population increases, presents a blue- 
shift as observed experimentally, We also show similar results with different either initial 
polarization (Fig. |^) or population (Fig. ^). The main difference between Fig. |] where 
P = 60% and Fig. |l| where P = 80% is the decrease of the ±1 splitting as well as a shorter 
duration of this splitting. In figure ^ we set the initial population na"^ equal to 0.05 (about 
5.0 10^° excitons per cm^) and the initial polarization once again P = 80%. The splitting is 
smaller than the one observed in the Fig. , about 0.25£'^^ and its duration is longer. 
The main consequences we extract from these results are the following: 

i) the spin level degeneracy breaking is proportional to the polarization of excitons pop- 
ulations and both polarization and splitting disappear in a time of the order of 50 ps 

ii) the ±2 excitons are negligibly populated and therefore do not influence the ±1 exciton 
energy levels although the opposite is true: the ±1 excitons do influence the ±2 energy levels. 

In order to compare with experimental information, an interesting way of presenting our 
results is to give energy levels at fixed times as a function of the initial exciton density. 
In figure ^ some of the conclusions we stated before become more clear. The splitting 
increases with P and n. We also observe how the red-shift of the spin —1 exciton at 
an early stage (lOps) becomes a blue-shift at 25ps. Those splittings have been observed 
experimentally but there is significant difference with our results. Some experiments 
|3^ seem to suggest that the majority carriers level remains essentially constant as a function 



of the initial density, something that clearly does not happens in our results. Moreover, our 
theory gives good results for level splittings but fails to describe the absolute position of 
the excitonic luminescence peaks as a function of the initial carrier density. We believe 
that these disagreements can be due to the pinning of excitons to impurities, an effect that 
might change the exciton energies and that has not been taken into account in our theory. 
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Another possible explanation for the mentioned disagreements could be the failure of the 
approximation of considering the center of mass at rest to explain non-resonant excitation 
experiments. 

In summary, we have presented a theory which describes the interacting polarized ex- 
citon gas with moderately high density. We have obtained a set of multi-exciton Wannier 
equations that we have solved, both perturbativately and variationally, obtaining blue-shifts 
originated by the inter-excitonic exchange interaction. Whenever differences between the ex- 
citon populations exist, energy levels are splitted. Calculated and experimentally observed 
splittings are in good qualitatively agreement. 
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APPENDIX A: CALCULATION OF h AND h 

In this appendix we calculate Ji and I2 both in 2D and 3D. The calculation of 




(Al) 



can be done both in 2D and 3D by using the Schrodinger equation 




(A2) 



Substituting ^2| in ^T] we obtain 




(A3) 



In the 2D case we obtain for this integral 




(A4) 
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and in the 3D case: 



CO t'^ rli' 

, , ^ =2l7r|igg^|a^ 



(A5) 



The calculation of 



dq 



df 



is more cumbersome. We start with the 2D case by perfoming the integral 



(A6) 



^00 



dt 



(2vr) 



ds 



ds 



(27r)2 



-y{s)e 



(A7) 



where 



dt 



1 f°° M , . m9 9 /""^ tJo(tr)dt 7^ 

2^i ^^o(tr)|0o(t)P^a^/^ ^ T 



^2(7) 



(A8) 



where Jo(x) is the Bessel function of the first kind, -^'2(7) is the Modified Bessel function 
and 7 = 2r/a. Since V{s) is the Fourier-transform the Coulomb potential, 



ds 



■V{s)e- 



(A9) 



and we obtain 



Eoo = ^ / r^K^l-) Jo(gr)t/r = e'a ^F^-, -, 1, -iqa/2f] (AlO) 
a'' Jo a 2 2 



where 



2Fi[a, 6, c, xj = 2^ —^^ — ; (0^) 



(a + fc- 1)! 
(a-l)\ 



(All) 



is the confluent Hypergeometrical function. In order to obtain I2 we use |A10| in and we 

get 



X 2Fi[\,l,l,-x^] dx 315 



[l + x 



2\3 



4096 



^3 I 77121? I 2 



(A12) 
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where we have used 



foo x2Fi[ll,l,-x^] dx _ 2 315 

Jo (l + x2)3 "3''4096 ^^^^^ 



(1 + X2)3 

In 3D the calculation of I2 is easier. We have 



(27r)'^'^"'^'' i (27r) 
= — e^a^ / —r / di / ^ J A14 

TT io (1 + X2)4y_i (1 + ^2 + - 2x|/0 

where x = a\q\, y = a|i| and a^q t — xycos{^). Now we perform the integrations over ^ and 
y in two steps 



r+l 

F{x,y) = 



d^ 



(1 + a;2 + 2/2 _ 2xy^y 
1 1 



6X1/ [ (1 + a;2 + _ 2a;|/)-^ (1 + a;2 + ^2 _^ 2x1/)' 

and 

dy ^. , TT 15 + 10a;2 + 3a;^ 



(A15) 



r°° dy ^, , TT 15 + lUx' + 6x" ^, , / , „x 



Hence, we obtain 



A 2''|f3D, 3 f^ xJ{x)dx 100 3^ 3 
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FIGURES 

FIG. 1. (a) ±1 exciton population as a function of time (0-400 ps). The dashed line is the 
— 1 population. In the inset we plot the +2 (circles) and -2 (crosses) populations (notice the 
different time scale), (b) Exciton levels as a function of time (0-140 ps) measured in quantum 
well free exciton energy (see text). Crosses and circles are again the ±2 excitons. The dashed line 
correspond to (-l)-exciton energy. Initial density na^ = 0.1, initial polarization P=80%. 



FIG. 2. As in fig. 1 with n=0.1, P = 60% 



FIG. 3. As in fig. 1 with n=0.05, P -- 

FIG. 4. ±1 exciton energies at (a) t = IQps and (b) t = 25ps as a function of initial densities. 
The solid and dashed lines correspond to an initial polarization of 80% and 60% respectively, and 
upper and lower lines correspond to spin +1 and spin —1 exciton respectively. 



24 



